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Abstract 



o 

5_l We introduce a tractable multi-currency model with stochastic volatility and correlated 

stochastic interest rates that takes into account the smile in the FX market and the evolution 
of yield curves. The pricing of vanilla options on FX rates can be performed efficiently through 
the FFT methodology thanks to the affinity of the model. Our framework is also able to 
describe many non trivial links between FX rates and interest rates: a second calibration 
exercise highlights the ability of the model to fit simultaneously FX implied volatilities while 
being coherent with interest rate products. 

1 Introduction 

C4-H The FX market is the largest and most liquid financial market in the world. The daily volume of 

^ FX option transaction in 2010 was about 207 billion USD, according to Mallo (2010). The stylized 

1— — 1 facts concerning FX options may be ascribed to two main cathegories: features of the underlying 

exchange rates and the implied volatilities respectively. The first and most important feature of 
FX rates is that the inverse of an FX rate is still an FX rate, so if S d ' f (t) is a model for EUR/USD 
exchange rate, i.e. the price in dollars of one euro, thus reflecting the point of view of an American 
investor, then S^' d (t) = 1/S d 'f(t) represents the USD /EUR rate, i.e. the price in euros of one 
dollar, hence representing the perspective of a European investor. This basic observation may 
be further generalized so as to construct e.g. triangles of currencies where then the no-arbitrage 

C\| relation S f ' d (t) — S^ k (t)S k ' d (t) must hold. This particular property of FX rates must then be 

coupled with the presence of a volatility smile for each FX rate involved in a currency triangle, in 
such a way that we can consider a model who is able to capture jointly relations among underlyings 
. . and their respective implied volatilities. 



Since the financial crisis, investors look for products with a long time horizon that are supposed 



to be less sensitive to short-term market fluctuations. Following Clark (2011), the risk involved in 
such structures may be intuitively understood in terms of a simple example. The value of an ATM 
call option, with maturity T, in a Bl ack-Scholes settin g, is usually approximated by practitioner s 



by means of the formula 0.04<7vT, see Bachelier (1900) and Schachermayer and Teichmann (2008). 
Under this approximation, the vega of such a position is simply given by 0.04V? 1 and hence scales 
as square root of T. If we look instead at the rho risk it can be shown that it scales as T. This 
heuristic observation suggests that when we look at shortdated FX options, the volatility smile 
is a dominant factor, whereas, as the maturity increases, the interest rate uncertainty plays an 
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increasingly important role. As a consequence, when we consider longdated FX products, a model 
that is simultaneously able to take into account volatility and interest rate risk is preferable. 



Besides the simplest case of currency swap, longdated FX products received increasing im- 
portance since the appearence of new structured products like power reverse dual currency notes 



(PRDC). According to Clark (20111, for some currency pairs "exchange rates in the future tend, 
once we actually get to those future dates, to be closer to current spot levels than the forward 
rates observed today would predict" . PRDC notes try to exploit this decoupling phenomenon so 
as to grant a speculative gain to the investor. PRDC notes pay coupons in case the future ex- 
change rates S d ^ {T k ) are closer to the current level of spot S d ''(t) rather than the forward price 
F{t,T k ) = S dJ (t)D d (t,T k ) / Di (t,T k ), where D l (t,T), I = d, f denote the discount factors under 
the currencies /, d respectively. Consider a sequence of payment dates {Tj , T/v}, where typically 
Tn is 30 years and AT k is six months. At these dates the following coupon is paid 

Vtu = AT fc max • 



where C^,C d are constants typically equal to respectively 0.15 and 0.1. The coupon may be 
rewritten as 

]\fd-f<f (~<d 

V Tk =AT kWlJ ^(S d 'f(T k )-K) + , K = S d >f{t)^j. (2) 

The basic example above can be further generalized so as to include capped or floored coupons. 
Futher examples of longdated FX products are callable, trigger or chooser PRDC notes. Given the 
long maturity of these products, it is essential to account for both oscillations in FX and interest 
rates. 



1.1 Related literature 



Historically, it has been quite a standard practice to employ, for the pricing of FX options, models 
originally designed e.g. for equity options, modulo minor changes. This approach was initiated 
by Garman and Kohlhagen] ( |1983 ), who adapted the Black and Scholes (1973) model, see Wystup 
(20061. Another important example is given by the Heston model, whose FX adaptation is pre- 
sented e.g. in|Clark| ( |2011[ ) and Janek et al. (2011 ). Other examples are the Stein and Stein Stein 



and Stein ( 1991 1 and the Hull and White Hull and White ( 19871 model, as described e.g. in Lipton 



(2001 1. Another recent contribution is Leung et al. (20131 where they adopted the Wishart multi- 



factor stochastic volatility model introduced in Da Fonseca et al 



setting. With a view towards quanto options Branger and Muck 
process. 



( 2008 ) in a single exchange rate 



(2012) also employ the Wishart 



The typical limitation of the previous approaches is that they neglect the relationships among 
multiple exchange rates. In fact, the joint presence of triangular relationships and volatility smiles 
makes it difficult to model all FX rates in a triangle of currencies. In the literature, there exists a 
stream of contributions that try to recover the risk neutral probability distribution of the cross ex- 
change rate, either by means of joint densities or copulas. These approaches represent an evolution 



of the Breeden and Litzenberger Breeden and Litzenberger ( 1978 1 approach. Among others, we 
recall Bliss and Panigirtzoglou (2002) and Schlogl (20121 and Austing (2011). Copulas have been 
employed in Bennett and Kennedy (20041 



Salmon and Schneiderl ( 2006 ) and Hurd et al 



In the presence of a stochastic volatility model of Heston type, Carr and Verma (2005) try to 



(2005). 



solve the joint valuation problem of FX options by specifying the dynamics of two rates influenced 
by a common stochastic volatility factor. This rather restrictive approach seems however difficult 
to extend. Asymptotic formulae for a SABR specification are provided in Shiraya and TakahasE] 
(2012 1. The approach we are interested in is presented in De Col et al. (2011 ), where a multifactor 



stochastic volatility model of Heston type is introduced. The model is coherent with respect to 
triangular relationships among currencies and allows for a simultaneous calibration of the volatility 
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surfaces of the FX rates involved in a triangle, like EUR/USD / JPY. The idea of De Col et al 
(2011 1 is inspired by the work of Heath and Platen (2006), who consider a model for FX rates of 
the form 



(3) 



where Q l , Qi represent the value of the growth optimal portfolio under the two currencies. A similar 
approach, known as intrinsic currency valuation framework, has been independently proposed in 



Doust ( 


2007) and 


Doust 


(2012) 



In the industry, the evaluation of longdated FX products is usually performed by coupling- 



Hull and White Hull and White ( 1993 ) models for the short rates in each monetary area with a 
stochastic process for the FX rate, which is usually assumed to be a geometric Brownian motion, 
see Clark (2011). Many authors, have studied the problem of combining uncertain FX rates or 
stocks with stochastic interest rates. Such a model represents a starting point for the approach 
presented in Piterbarg ( 2006[ ), where a local volatility effect in the FX process is also introduced. 



Van Haastrecht et al. (2009) derive closed-form pricing formulae under the same kind of three 



factor model, coupled also with a stochastic volatility process of Schobel and Zhu ( 1999 ) type, see 



also Ahlip (2008) and Ahlip and King (20101. This kind of setting implies that both interest rates 



and volatility may become negative. Van Haastrecht and Pelsser (2011) extend the previous ap- 



proach and also consider approximate solutions when the instantaneous variance follows a square 



root process. Deelstra and Rayee (2012) first present a local volatility framework, and then extend 



it to a stochastic volatility setting and provide approximations by means of a Markovian projection. 

A desirable feature, when we couple interest rates and FX rates with stochastic volatility, is 
that we would like to use models such that the instantaneous variance and the interest rates remain 
positive. A natural choice in this sense is given by the square root process. However, introducing 
a non-zero correlation among square root processes breaks the analytical tractability of the model, 
when we work with affine processes on the canonical state space M> x R n . IGrzelak and Oosterlee 



(2011 1 attack the problem by providing approximations of the non affine terms in the Kolmogorov 



PDE satisfied by the characteristic function of the forward exchange rate. 



1.2 Main results of the article 



In this paper we propose an extension of De Col et al. (2011) to the case where the stochastic 



factors driving the volatilities of the exchange rates belong to the matrix state space Sj", the cone 
of positive semidefinite d x d matrices. Stochastically continuous Markov processes on with 
exponential affine dependence on the initial state space have been characterized in Cuchiero et al 
( 2011[ ), see also Mayerhofer (2012 1. We will see that all good analytical properties of the De Col 



et al. (20111 approach are preserved, while we allow for a more general dynamics. Our model 



is at the same time an affine multifactor stochastic volatility model for the FX rate, where the 



instantaneous variance is driven by a Wishart process, see Bru (1991), Da Fonseca et al. (2008), 



(20081, Buraschi et al. (2008), Chiarella et al. (2010) and Gnoatto (20121. The model has many 



and a Wishart affine short rate model, see Gourieroux and Sufana| ( 2003[ ), Grasselli and Tebaldi 



interesting features, namely: 

• it can be jointly calibrated on different FX volatility smiles; 

• it is coherent with triangular relationships among FX rates; 

• it allows for closed form solutions for both FX options and basic interest rates derivatives; 

• it allows for non trivial correlations betwen interest rates and volatilities. 



The first two interesting features are shared with the model in De Col et al. (2011). The other 
results constitute a novel contribution of the present article, which is outlined as follows: In 
section [2] we set up our modelling framework and provide an example highlighting the flexibility 
of the approach. Section [3] is devoted to a complete characterization of all risk-neutral measures 
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associated to the different economies whiie Section [4] shows that the present setting allows for 
stochastic correlations among many economic quantities. Section [5] presents all closed form pricing 
formulas for FX and interest rates derivatives, together with asymptitic expansions of FX implied 
volatilities. The numerical treatment of our model is the topic of Section |5j where fit the model to 
an FX surface and two yield curves. Finally Section [8] summarizes our findings. Technical proofs 
are gathered in the Appendix. 



2 The model 

We consider a foreign exchange market in which N currencies are traded between each other 
via standard FX spot and FX vanilla option transactions. The value of each of these currencies 
i = 1,..,N in units of a universal numeraire is denoted by S°' l (t) (note that S°' l (t) can itself be 
thought as an exchange rate, between the currency i and the artificial currency i = 0: we will see 
that the results are independent of the specification of the universal numi'^raire). 

We assume the existence of N money-market accounts (one for each monetary area), whose 
values are driven by locally deterministic ODE's of the type: 



dB\t) =r l (t)B l (t)dt, 



i = l, 



,N: 



(4) 



and we denote with r° the interest rate corresponding to the artificial currency. 

We assume the existence of a common multivariate stochastic factor S driving both the interest 
rates and the volatilities of the exchange rates S ' 1, (t) . We model the stochastic factor S as a matrix 
Wishart process (see Bru ( |1991 )) evolving as 



dS(t) = (flfl T + ME(t) + Y,(t)M T )dt + y/E(i)dW(t)Q + Q T dW(t) T . 



(5) 



We assume ft, M, Q e Md, the set of d x d real matrices, W = (Wt) t>0 € Md is a matrix Brownian 
motion (i.e. a dx d matrix whose components are independent Brownian motions). The dimension 
d can be chosen according to the specific problem and may reflect a PCA-type analysis. In order 
to ensure the typical mean reverting behavior of the process we assume that M is negative semi- 
definite, moreover the matrix 12 satisfies the condition f251 T = f3Q T Q, for j3 > d + 1. This last 
condition ensures the existence of a unique strong solution to the SDE ([5| , according to Corollary 
3.2 in |Mayerhofer et al.| ( [20Tl] ). 

We model each of theS^i) via a multifactor Wishart stochastic volatility model. Formally, 
we write 



dS°' l (t) 



(r°(t) - r\t))dt - Tr \A % y/Yfi)dZ{t) 



1. 



(6) 



where Z t € Md denotes a matrix Brownian motion. 

The diffusion term exhibits a structure that is completely analogous to the one introduced in 



Heath and Platen (2006) and De Col et al. (2011 ): in the present case we have that the dynamics 



of the exchange rate is driven by a linear projection of the variance factor -^/S(t) along a direction 
parametrized by the symmetric matrix Ai. As a consequence the total instantaneous variance is 
Tr [AiT,(t)Ai] dt. 

The stochastic factor S drives also the short interest rates: 



h° + Tr [J?°E(t)] 
h l +Tr [H^{t)] , 



(7) 

(8) 



for h k > 0, H k eS+,k = 0,..,N. 

We assume a correlation structure between the two matrix Brownian motions Z{t) and W(t), 
by means of an invertible matrix R according to the following relationship: 



W{t) =Z{t)R v + B(t)y/ Id — RR T , 



(9) 
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where B(t) is a matrix Brownian motion independent of Z(i). We denote by S l,3 (t), i,j = 1, .., N 
the exchange rate between currency i and j. By Ito's lemma we have that S % ' 3 (t) = S ' 3 (t)/ S°' l (t) 
has the following dynamics 



dS l > 3 (t) 



= (r\t) - r j (t))dt + Tr[{A % - Aj)E(t)Ai]dt 
+ Tr[(A l -A J )^(t)dZ(t)}. 



(10) 



The additional drift term in ( 10 ) can be understood as a quanto adjustment between the currencies 
i and j. 

We try to provide some intuition concerning the flexibility of the approach by considering an 
introductory example. 

Example 1. Let d = 4. Consider the case 



= h l + Tr 



r 3 = hp + Tr 
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E(t) 



E(i) 



= h l + Z n (t) 



= h j +X 2 2(t) 



Moreover we let the matrices A l ,A J be partitioned as follows 



A' 



/ \ 



A\ x A\ 2 

V A\ 2 A 22 ) 



The idea is that the first two diagonal elements of the Wishart process are mainly responsible 
for the dynamics of the short-rates whereas the third and the fourth drive the volatilities. The 
elements of the Wishart process are related among each other in a non-trivial way, in fact 



d<E ii ,E ifc )t = (s« (Q T Q) Jt + £<i (Q T Q) Jk + Eji (Q T Q)i 



dt. 



(11) 



Consequently, in the present four- dimensional hybrid FX-short-rate model, the correlation between 
e.g. the element En (which is mainly responsible for the dynamics of the first short rate) and E 33 
( which drives the volatility of the FX rate ) is stochastic and given by 

d(E 11 ,E 3 3) t = 4E 1 3 (Q T Q) 13 dt. 

From the discussion above, we realize that we are in presence of a hybrid model, allowing for 
non trivial links among interest rates and FX rates, while preserving full analytical tractability, as 
we will see in the sequel. This interesting feature, to the best of our knowledge, is not shared by 
other existing hybrid models, that usually rely on approximations of the solution of the Kolmogorov 
PDE involving non-affine terms (see e.g. Grzelak and Oosterlee (2011), Grzelak et al. (2012}) ). 



3 Risk neutral probability measures 

Up to now we have worked in the risk neutral measure defined by our artificial currency. In practical 
pricing applications, it is more convenient to change the numeraire to any of the currencies included 
in our FX multi-dimensional system. Without loss of generality, let us consider the risk neutral 
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measure defined by the i-th money market account B l (t) and derive the dynamical equations for 
the standard FX rate S l,: >(t), its inverse S J ' l (t), and a generic cross S J ' l (t). 

The Girsanov change of measure that transfers to the Q l risk neutral measure (i.e. the risk 
neutral measure in the z-th country) is simply determined by assuming that under Q 1 the drift of 
the exchange rate S % ' J (t) is given by r l — r 3 (or equivalently by the fact that the money market 
account B 3 (t) is a Q'-martingale once discounted by the interest rate r*(i)). The associated 
Radon-Nikodym derivative is given by 



: exp 



Tr[A % ^(s)dZ{s)] 



Tr[A i T,{s)A i }ds 



In Mayerhofer ( 2012 ), conditions ensuring that the stochastic exponential above is a true martingale 
are provided. In the following, we proceed along the lines of De Col et al. (2011 ). The possibility 
of buying the foreign currency and investing it at the foreign short rate of interest, is equivalent 
to the possibility of investing in a domestic asset with price process B* = B 3 S 1 ' 3 , where i is the 
domestic economy and j is the foreign one. Then 

dB){t) =d(B 1 (t)S l ' 1 (t)) 

= B 3 {t)S l ' 3 {t) {{r l {t) - r 3 {t))dt + Tr[(A t - A^^A^dt 

+Tr[(A l - A J )y r JMf)dZ(t)}^j + li'S' ■ihr'Uull 

= B){t) (r\t)dt + Tr[(Ai - A j )Y.(t)A i \dt + 7V[(i4, - A^y^dZit)} 
)(t) (r\t)dt + Tr[{A t - A i )^/S(^)^ Q *(^)]) , 



Bl 



where the matrix brownian motion under Q l is given by 

dZ® =dZ + y/^t)Aidt, 

then the Q'-risk neutral dynamics of the exchange rate is of the form 



d5* J (t) 



B 3 {t) 



S h3 (t) ((r'(t) - r 3 {t))dt + Tr [(A, - A,) 03(4)^ (t)] ) 



The measure change has however also an impact on the variance processes, via the correlation 
matrix R introduced in The component of dB{t) that is orthogonal to the spot driver dZ(t) 
is not affected by the measure change; this is a natural choice that is consistent with the foreign- 
domestic symmetrj|^| We are now able to derive the risk neutral dynamics of the factor process 
£(i) governing the volatility of the exchange rates under Q 1 , that is given by 



(12) 



dW q \t) = (dZ(t) + VZAidtj R T + dB{t)y/l d - RRJ . 

From ([5| and ^ we derive the Q l -risk neutral dynamics of £ as follows: 

dT,(t) =(fifi T + ME(t) + Y,{t)M T )dt 

+ ^/nj)((dZ{t) + ^/ntjAidt) R T + dB(t)y/l d - RR T ) Q 



Q T \^R (dZj + A l ^/jMt)dt S j + yjl d - RR t1 dB T J y^(f) 
Y,(t)A i R T Qdt ~ Q T RA l Z{t)dt. 



1 This assumption is in line with the procedure that has be en introduced in|De Col et al.| ( |201lfr and implies that 
the model is consistent with the foreign-domestic parity as in |Del Baiio Rollinj ( |2008^ 



G 



Now define 



M -Q RAi 



so that using ( 12 ) we can finally write 

d£(t) = {nn T + M Q *Y,(t) + Y,(t)M Q *< T )dt 

^W)dW Q \t)Q + Q T dW^- T (t)^({), 
from which we deduce the relations among the parameters: 

= R, 

Q Qi = Q, 

= M-Q T RAi. 



(13) 



We observe that, like in the multi-Heston case of De Col et al. (2011), the functional form of 
the model is invariant under the measure change between Q u and the ith-risk neutral measure. 
The inverse FX rate under the Q'-risk neutral measure follows from Ito calculus, recalling that 
Si,i = IS' -' ) ': 



dS^jt) 



= S l > 3 (t)d 



= (f(t) - r*(t) + Tr [{A, - A^t^A, - AA}) dt 



Tr 



(A, - M)\/W)dZP{t) 



which includes the self-quanto adjustment. Similarly, the SDE of a generic cross FX rate becomes 
dS*< l {t) _S l >i{t) fS^it) 



(r j (t) - r\t) + Tr [(A, - 4)E(t)(A,- - At)]) dt 



Tr 



The additional drift term is the quanto adjustment as described by the current model choice. 
By applying Girsanov's theorem again, this time switching to the Q J risk neutral measure, the 
term is removed while the Wishart parameters change according to the following fundamental 
transformation rules: 



Q Q3 =Q q \ (14) 
M QJ = Af Ql - <5 Qi,T i? Qi (Aj - A{). 

4 Features of the model 

4.1 Functional symmetry of the model 

Recall that a crucial property of the FX market requires that products or ratios of exchange rates 
in a triangle are also exchange rates, meaning that the dynamics of the exchange rates must be 
functional symmetric with respect to which FX pairs we choose to be the main ones and which 
one the cross. That is, it is not a priori trivial to obtain a model such that the dynamics for the 
inverse of an exchange rate shares the same functional form in its coefficients. Or equivalently, the 
dynamics of (S^S 1 **) computed by applying the Ito's rule to the product (S^/S ' 1 ) x (S°' j /S ' 1 ), 
must give the dynamics of a process that shares the same functional form of both S 1 ' 1 and S 1 ^ . 
This symmetry property is fundamental in order to be able to joint calibrate and price consistently 



multi currency options (see e.g. De Col et al. (2011 1). 
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Proposition 1. The dynamics of the exchange rates satisfies the triangular relation, namely 
the model is functional symmetric. 



Proof. Sec the Appendix. 
4.2 Stochastic Skew 



□ 



In analogy with De Col et al. (20111, if we calculate the infinitesimal correlation between the log 
returns of S 1,3 and their variance Vol 2 (S 1,3 ), we find that it is stochastic. This is a nice feature of 
the model since it implies that the skewness of vanilla options on S l,J is stochastic, which is a well 
known stylized fact in the FX market (see e.g. Carr and Wu (2007)). 
Let us consider the infinitesimal variance of S 1 ' 3 : 



d(S i ' j ,S i ' j ) t = Tr[(Ai - Aj)X(t)(Ai - Aj)]dt, 



so that we can write 

SW(t) 

where we defined 



= (r l - r j )dt + Tr \{A l - Aj ) y^YjtjdZ®' (t) 



(r 4 - r>)dt + ^Tr [(A, - A^t^A, - A j )]dB 1 (t), 



(15) 



dB x (t) 



Tr 



(A i -A j )^)dZ^(t) 



^/Tr[(A l -A J )m(A-A J )] 



The process B\ — (Bi(t)) t>0 is a scalar local martingale with quadratic variation given by t, hence 
an application of Levy characterization theorem allows us to claim that Bi is a Brownian motion 
and we still denote by S 1,3 the weak solution to the SDE (15) . 
The dynamics of the variance is given by: 

dTr [(Ai - Aj)T,(t)(Ai - Aj)] 

= (Tr [{A, - A J )nn T {A l - Aj)] + 2Tr [(A l - Aj)MT,(t)(Ai - Aj)]) dt 
+ 2Tr \(Ai - Aj)\/Y,dW Qi {t)Q{Ai - Aj) . 

In order to determine the scalar Brownian motion driving the variance process we shall compute 
the following quadratic variation: 

d(Tr[(A i -Aj)^(A i -A j )]) t 

= 4d Q Tr [(A, - Aj)^Jnv)dW® ' (v)Q(A t - Aj)] , 
Tr \(Ai - A j )y^E(ujdW Qi (u)Q(Ai - Aj) 

d 

£ (Ai - Aj) ah jnv) bc dWZ (v)Q de (Ai - Aj) e 

a,fe,c,d,e— 1 







4d< 



(A - ^^ysM^^WQ^A, - Aj)i p 

■,q,r,s,Z— 1 / 
d 

(Ai- A j ) ab yjYfi) br Q se (A i - Aj) ea 

x (Ai - Aj) pq y/Wf) qr Q s i(Ai - Aj)i p dt 
= ATr [{A, - Aj) 2 E(t)(Ai - A j ) 2 Q T Q] dt. 
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Then we can use the same arguments as before and express the dynamics of the variance as follows: 
dTr [(Ai - A j )'E(t)(A i - Aj)] 



(...) dt + 2JTr [(A, ~ AtfTftXAi ~ Aj) 2 Q T Q]dB 2 (t), 



where 



dB 2 (t) := 



Tr [(A t - AjyVXdW® \t)Q{Ai - Aj) 
^Tr [(A, - - A 3 )WQ] 



which allows us to compute the covariation between the two noises. Notice that in the calculation 
above we are assuming the invariance of the correlation with respect to the change of measure that 
was explained in Section [3] The skewness is then related to the quadratic covariation between the 
two noises Bi,B 2 : 



d(B 1 ,B 2 ) t 



d (j- Tr \(Ai - A^^nfidZf] , J Q Tr \(A Z - A^y/^jdW^ (t)Q(Ai - A,) 



^Tr [{A, - A 3 )^(t)(A t - Aj)]y/Tr [(A, - A 3 ) 2 Z(t)(A 
d (/ Tr [(A, - A^jWfjdZ® {t)\ , f Q Tr [(A t - A 3 )^Z(tjdZ® (t)R T Q(A, - A,) 



^jTr [(Ai - AjMtHAi - Aj)]y/Tr [(A t - A 3 ) 2 X(t)(A* - A) 2 Q T Q] 

Y.'!,.,,,...,.. : (A ~ A^yJWti^dZftWAi - A 3 ) ab ^W) b dZ%(t)R T de Q ef (A^ - A 3 ) fa 

^jTr [(Ai - Aj)T,(t)(Ai - A 3 )]^Tr [(A, - A 3 ) 2 Z(t)(A t - A 3 ) 2 Q T Q] 

>::;,/,.,,./,, : a - a^j^^j^j^ - a^ua - A 3 ) af Qj e R ed 

^Tr [(A, - Aj)^(t)(Ai - Aj)]y/Tr [(Ai - A^WMi - A 3 ) 2 Q T Q] 



dt 



Tr [(Ai - A 3 )E(t)(A, - Aj) 2 Q T R] 



^Tr [(A, - Aj)E(t)(Ai - A 3 )]^/Tr [(Ai - A^WM^i - A 3 ) 2 Q T Q 



dt. 



(16) 



This quantity is proportional to the skew: in particular, by looking at the numerator we realize 
the proportionality in the asymptotic expansion of the Proposition [5] 

4.3 A stochastic variance-covariance matrix 

We would like to discuss the positive defmitcness of the variance-covariance matrix. For simplicity, 
we consider the case of three currencies, meaning that we will have a 2 x 2 candidate covariance 
matrix: 



We know that 



(\nS i ' i AnS i ' l ) t (lnS" M ) t 



d(\nS^) t = Tr [(A z - Aj) £(*) (Ai - A,)} dt, 
d (In S i ' j ,lnS i ' l ) t = Tr [(Ai - Aj) E(t) (A, - At)] dt. 



(17) 



(18) 
(19) 



We first look at (18 1. We recall that we assumed Ai,Aj,Ai £ S d - Without loss of generality 
(otherwise put V = -V for V £ Sj), let (A, t - Aj) £ Recall that the cone Sj is self dual, 
meaning that: 

S+ = {ueS d | Tr[uv}>0,yveS+}. 

Let O be an orthogonal matrix, then we may write: (Ai — Aj) = OAO T , where A is a diagonal 
matrix containing the eigenvalues of (Ai — Aj) on the main diagonal. Then we have: 

Tr [(Ai - Aj) E(t) (Ai - Aj)] = Tr [OKO T Y J (t)OW T ] 

= Tr [£(i)0A 2 O T ] > 
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by self-duality. This shows that variances are positive. Now we would like to check that the 

. Now let 



variance-covariance matrix is in Sf . Now let 



AT(t) = (Ai - Ai) v/E(t), 

then, using the Cauchy-Schwarz inequality for matrices we have 

Tr [(Ai - A,) E(t) (Ai - A,)} Tr [(A, - A t ) E(t) (A, - A l )\ 
= Tr [M(t)M T (t)] Tr [Af (t)J\f T (t)] 



> 



Tr[M(t)N T (t)} 
= Tr[(A i -A j )V(t) (Ai ~ Ai)] 2 . 

This implies that the determinant of the instantaneous variance-covariance matrix is positive, so 
we conclude that the variance-covariance matrix is well defined, and, as a side effect, we have the 
usual bound for the correlations, i.e. all correlations are bounded by one (in absolute value). 

4.4 Stochastic correlation between short rates and FX rates 

We now show that our model allows for non trivial dependencies between the short rates and the 
exchange rates. 

Proposition 2. The instantaneous correlation between a short rate r % and the log-FX rate log S lJ (t) 
is stochastic and given by 

= Tr [(A l -A J )Z(t)HQ T R\ m 

/VjOSSI ' 3( * ) V Tr Wi - A,) kt) (Ai - A,)WTr p*Ep5gT] " 

Proof. Under the Q l -risk neutral measure, the covariation between the short rate and the log-FX 
rate is given by 

d (log S lJ ,r l ) t =dl^J^Tr [(^ - Aj) y/Z(u)dZ(u) , jf 2Tr ^Hy/i:(u)dZ(u)R T Q J 

= 2d( ^ / (Ai- Aj) i:j o-j k (u)dZf.i(u), ^ H pq a qr (u)dZ rs (u)R ns Q np J 

\ij,k=l^° p,q,r,s,n=l l t 

— (Ai Aj) 3 j G j r G r qH qpQ n pR ns dt 

= Tr [(Ai - Aj) T,(t)HQ T R] dt, 
while for the short rate we have 

d (r\ r l ) t = Tr [QWZ(t)H l Q T ] dt. (21) 

Given the instantaneous quadratic variation of the log-FX rate and the short rate we conclude. 

□ 

4.5 Stochastic correlation between short rates and the variance of the 
FX rates 

The richness of our model specification may be further appreciated when we look at the correlation 
between any of the short rates and the variance of the FX rates, which is not usually captured in 



the literature (see Grzelak et al. ( 2012 1 ) . 
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Proposition 3. The instantaneous correlation between the short rate r % and the variance of the 
FX rate is stochastic and given by 



Tr [(A, - A,) T.(t)H i Q T Q (A, - A,)] 



Tr 



(22) 



np 



Pr z .Var l >3 

(A, - Aj) 2 E(t) (Ai - Aj) 2 g T Qj y/Tr [QH i Y.(t)H i Q T ] 

Proof. Recall that the noise of the scalar instantaneous variance process is 

dTr \(A l - Aj) £(i) {A, - Aj)] = (...) dt + 2Tr [(A, - Aj) ^/WfijdWQ {A, - Aj 
and observe that the noise of the short rate is simply given by 

dr l (t) = (...) dt + 2Tr liP y/ltydW(t)Q 

consequently 

d{r\Tr[{A i -Aj)H{A i -Aj)]) t 

= d \J 2Tr \H i y/T,{u)dW{u)Q , 2Tr [(A, - Aj) y/H(u)dW(u)Q {A, - A 

I „ d r. , d 

= d \ X! H ij°~jk(u)dWki(u)Qu, {Ai- Aj) pq a qr (u)dW rs (u)Q sn (Ai- Aj) 

V° i,j,k,l=l p,q,r,8,n=l I t 

= 4 HijO-j r (t)Q si (A, - Aj) m a qr (t)dW rs (t)Q sn (Ai - Aj) tp dt 

i,j,p,q,r,s,n=l 

= 4 J2 A i) Pq ^r{^rj{t)HjiQj s Q sn (Ai - Aj) np dt 

i,j,p,q,r,s,n=l 

= Tr [(Ai - Aj) Y,(t)HQ T Q (A, - Aj)] dt. 

Combining the above results with the quadratic variation of the instantaneous variance and 
the quadratic variation of the short rate we conclude. □ 

5 Pricing of derivatives 

A distinctive feature of our model is the ability to price in closed form derivatives written on differ- 
ent underlyings, meaning that we can consider, in a unified approach, more markets simultaneously. 
In particular, we can consider jointly the FX and the fixed-income market. Basic European prod- 
ucts like calls on FX rates and interest rate products related to different monetary areas can be 
analyzed together in a single model. In principle, this feature allows the desk to be jointly fitted 
to different markets by means of a single model. In view of this, we provide fast pricing formulas 
in semi-closed form, up to Fourier integrals. 

5.1 European FX options 

We first provide the calculation of the discounted Laplace transform and the characteristic function 
of x l 'i(t) := lnS' JJ (i), that will be useful for option pricing purposes. Let us consider a call option 
C(S i ' j (t),K i ^,r),i,j = l,..,N,i ^ j, on a generic FX rate S^(t) = cxp(x 1 ^ (t)) with strike , 
maturity T (t = T — t is the time to maturity) and face equal to one unit of the foreign currency. 
For ease of notation set: B9' = R and Q Ql = Q and we will use the shorthand M Q * = M. 
We proceed to prove the following: Being an affinc model, the discounted characteristic function 
conditioned on the initial values 



:j (w,t,r,x,s) =Ef [ e -# r ' d V u * i ''( T V'- 7 '(*) =*,£(*) =£] 



(23) 
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can be derived analytically (here i = \J— 1). Standard numerical integration methods can then 
be used to invert the Fourier transform to obtain the probability density at T or the vanilla price 
via integration against the payoff, with overall little computational effort. In fact, from the usual 
risk-neutral argument, the initial price of the call option can be written as (domestic) risk neutral 
expected value: 

e -f t T rlds ^(T) 



C(S i > 3 (t),K i < 3 \r)=E^ 



and by applying standard arguments (see e.g. |Carr and Madan ( 1999 ) and also Bakshi and Madan 



(20001, Dufhe et al. (2000) and Sepp (2003)) it can be expressed in terms of the integral of the 
product of the Fourier transform of the payoff and the discounted characteristic function of the 
log-asset price: 



where 



1 



2tt 

e ±\ x r e x 



(f> i ' j (-i\,t,T,x,J:)<i>(X)dX, 



(24) 



K iJ ) + dx 



$(A) = 

is the Fourier transform of the payoff function and Z denotes the strip of regularity of the payoff, 
that is the admissible domain where the integral in ( 24 ) is well defined. In other words, the pricing 
problem is essentially solved once the (conditional) discounted characteristic function of the log- 
exchange rate is known. We recall the relationship between the characteristic function and the 
moment generating function. If we denote via G 1 ' 3 (uj, t, r, x, E) the discounted moment generating 
function, given by 



G lJ (w,i,T,x,E) =Ef [e-^ Tr = ds e^ ,3(T) |x JJ (i) = x,E(t) = E], 



(25) 



we simply have <p 1 ' 3 '(uj, t, r, x, V) = G 1 ' 3 (iw, t, r, x, E). As a consequence, in our affine model it 
is sufficient to derive the discounted Laplace transform, which is explicitly given by the following 
proposition. 

Proposition 4. Assume that the matrices Ai,Aj are symmetric. Then, in the Wishart model, the 
discounted Laplace Transform ofx t ' 3 {t) := log S l ' 3 (t) is given by: 



where: 

A= (uj (h l - hP) -h l )r- -Tr 

B(t)=B 22 (t)- 1 B 21 (t) 

and £>22(t), 621(1") are submatrices in: 

' Su(t) B 12 (t) \ _ 
Bai(rj B 22 (t) J 

M + ujQ T R(Ai - Aj) 



G i,j (u, t, T, x, E) = exp [ux + A(t) + Tr [B(r)E]] 

13 r 



log J-(r) + (il/ T + lo {A, - Aj) R T Q^j 



(26) 

(27) 
(28) 



expr 



(Ai - A } y + (w - l)W - loH 3 



M 1 



-2Q T Q 

- uj (Ai — Aj) R T Q 



Proof. See the Appendix. 



(29) 



□ 



As we stated above, for uj — iA (where i = \f—l), we obtain the discounted characteristic 

ate, hence 

exp {—ac} 



function of the log-exchange rate, hence we can compute call option prices as e.g. in Carr and 

+00 



Madan ( 1999 ) via: 



2tt 



$l{e- ivc tp(v)}dv, 



where: 



tp(v) 



G(±(v-(a + l)i),t,T,x,V) 
(a + iv) (1 + a + iv) 



(30) 



(31) 
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5.2 Expansions 



The calibration of our model relies on a standard non-linear least squares procedure. This will 
be employed to minimize the distance between model and market implied volatilities. The model 
implied volatilities are extracted from the prices produced by the FFT routine. This procedure 
for the Wishart model is more demanding from a numerical point of view than the analog for the 
multi-Heston case e.g. in De Col et al. (2011 ). An alternative approach is to fit implied volatilities 
via a simpler function. A possibility is to find a relationship between the prices produced by 
the model, and the standard Black-Scholes formula. The next result states that it is possible 
to approximate the prices of options under the Wishart model, via a suitable expansion of the 
standard Black-Scholes formula and its derivatives, analogously to what has been done in |De CoI| 



et al. (2011 1. The proof, that is reported in the appendix, relies on arguments that may be found 



Lewis (2000) and Da Fonseca and Grasselli ( |2011| (we drop all currency indices, it is intended 



that we are considering the (i,j) FX pair). Define r := T — t and let us define the real deterministic 
functions B°, B\ B 20 , B 21 as follows: 



B° = 


/V- 




Jo 






B 1 = 






Jo 




+ {Ai- 






B 20 = 


/V- 




Jo 






B 21 = 






Jo 




+ (Ai- 



e (r-«)MT ^. _ A .j e {r-u)M d ^ 



B°(u)Q T R(A l - Aj 



)M 



du, 



(b 1 {u)Q t R(A 1 -A, j ) 



Moreover, the real deterministic scalar functions A (t) , A 1 (t) , A 20 (t) , A 21 (t) are given by: 



A°(r)=Tr 
A 1 (r)=Tr 
A 20 (T)=Tr 
A 21 {r)=Tr 



an 1 
nn T 
nn T 



B°(u)du 
B x {u)du 
B 20 (u)du 
B 21 {u)du 



Finally, let 



be the integrated variance. 



a 2 r 



A°{T)+Tr S°(t)E 



(32) 

(33) 
(34) 

(35) 

(36) 
(37) 
(38) 
(39) 

(40) 



Proposition 5. Assume that all interest rates are constant and the vol of vol matrix Q has been 
scaled by the factor a > 0. Then the call price C(S(t), K, t) in the Wishart-based exchange model 
can be approximated in terms of the vol of vol scale factor a by differentiating the Black Scholes 
formula CbszS (S(t), K,a,r) with respect to the log exchange rate x(t) — In S(t) and the integrated 
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variance v = o t: 

C(S(t),K,r) «C B&S (S(t),K,a,r) 



B x (r)E(t) )d 2 xv C BSz s(S(t),K,a,r) 



Tr 



+ a z (A M (r)+Tr B zu (r)Z(t) )d^C BSlS (S(t), K, a, r) 



a 2 (A 21 (t) + Tr [Z? 21 (r)£(i)] ) S^ v C B k3 (S{t),K, a, r) 
(i 1 (r) + Tr [B^rjEft)] ) * ^C BkS {S{t),K, a, r) 



a 

T 



Proof. See the Appendix. 



(41) 

□ 



Finally, we can state another formula, that does not involve the computation of option prices 
and constitutes an approximation of the implied volatility surface for a short time to maturity. 
This formula may be a useful alternative in order to get a quicker calibration for short maturities. 

Proposition 6. Assume that all interest rates are constant. For a short time to maturity the 
implied volatility expansion in terms of the vol-of-vol scale factor a in the Wishart-based exchange 
model is given by: 



tTr [{Ai - Aj) E(t) {A i -A j )]+, 



Tr 



(A l -A J ) 2 Q T R(A l -A J )^(t) 
TrUi-Aj)^) (A, -A,)} 



m f 



Tr[(Ai-Aj)E(t) (Ai -Aj)Y 



-Tr [(A, - Aj) Q T Q (A, - A,) S(t)] 



-Tr 



[(A, - A,) Q T R (Ai - A,) + (Ai - Aj) R T Q (A, - A,)] 



x Q T R(A i -A j )E(t) 



5 Tr [(A, - A,) Q T R(A t - A,) E(t)]* 
4" 



where mf = log I — L -^ 

Proof. Sec the Appendix. 



Tr\(Ai- A^t) (Ai — Aj)] 
denotes the log-moneyness. 



(42) 



□ 



5.3 The pricing of zero-coupon Bonds 

The flexibility of the Wishart hybrid model opens up the possibility to perform a simultaneous 
calibration of interest rate and foreign exchange related products. This is a preliminary step 
that guarantees a coherent framework for the evaluation of payoffs depending on several interest 



Mercurio 


(2006 


) and 


Clark ( 


2011) 



bond, that constitutes a building block for many other linear interest rate products. The following 
proposition may be easily proved along the same lines of Proposition [4] 

Proposition 7. The price of a zero-coupon bond at time t, with maturity T , under a generic risk 
neutral measure Qd, *s given by 



,-f t T (h d +Tr[H d S(s)})ds 



exp {A zc (t) + Tr [B zc (r)E(i)] } 



where the deterministic functions A ZC ,B ZC satisfy the system of matrix ODE 
i>AZ( --Tr [PQ t QB zc (t)} - h d , A zc (0) = 0, 
--B ZC (t)M®« +M q *- t B zc (t) 



8t 
dB zc 
dr 



2B zc (t)Q t QB zc (t) - H d , B zc (0) = 0, 



(43) 



(44) 



(45) 
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whose solution is given by 

A zc (r) - ' t- fi.,,. r> z > ' 

and B22 1 (j\B>2i (t) are submatrices in: 



^Tr [log B| 2 C (r)+r 



exp 



-2Q 1 Q 



(46) 
(47) 

(48) 



Given this simple formula for zero-coupon bond, we may consider a joint calibration to the FX 
smile and to the risk-free curve of two different economies, so as to capture simultaneously the 
information coming from the FX smile and the interest rate curve, that plays an important role for 



long maturities. More precisely, a direct inspection of formula (43), reveals that the Wishart short 



rate model belongs to the class of affine term structure models, meaning that the yield curve is of 
a particularly simple form. For fixed we define the yield curve as the function Y : M>o — > R>o 
with 



Y(t) =-- (A ZC (T)+Tr [B zc '(r)S(t)]) 



(49) 



In Section 7. 1 we provide a concrete calibration example involving the simultaneous fit of the FX 
surface and two yield-curves. 



5.4 The pricing of a Cap 

For the sake of completeness, we also report the pricing of an interest rate cap. We would like 
to point out that the formula presented here is derived under the assumption of a single curve 
interest-rate framework. A generalization involving the recent developments in the context of in- 
terest rate modelling, like multiple-curve models/OIS discounting, is beyond the scope of the paper. 

Suppose we have a sequence of resetting dates {T 7 , • • • , Tp_i} and payment dates {T 7+1 , ■ ■ • , T«} 
at which the Libor rate L(Ti, Ti, T,+i) = L(T, Tj+i) observed at time T for the time span between 
Ti and Tj+i is paid. The price at time t < T 7 of an interest rate cap with notional N, strike price 
K is given by 



C(t,K)= NT * £(id [e-'« Tir « ds i;) -tfT 

i=7+l 



(50) 



It is well known that the expectation above may be conveniently rewritten as a put option on a 
zero coupon bond 



Cap(t,K)= N(l+nK) 

i=7+l 



x E K 



■ft 



1 + T,K 



7 -P(T i _ 1 ,T l ) 



(51) 



Given the nice analytical properties of the Wishart process, we would like to compute this 
expectation using Fourier methods. The standard approach, in the presence of such an expectation, 
would involve a change to a forward risk-neutral measure. While this is also possible in this case, we 
remark that the resulting dynamics of the Wishart process would involve time varying coefficients, 
which in turn would lead to the solution of time dependent Matrix Riccati differential equations, a 
highly non trivial problem. Fortunately, this kind of problem can be conveniently overcome, since, 
in our setting, it is easy to compute expectations that involve the discount factor and the zero- 
coupon bond. In fact, without loss of generality, let Ti = r, Vi and introduce a positive constant 
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a > and k = log( 1 2 rK )'- from e.g. Carr and Madan (19991 it follows that the price of a put on 
a zer-coupon bond can be written as 



p ak poo 

Put a (t,k) = —j R{e- ivk ^(v)}dv 



where 



ip{v) = 



-Jt Ti_1 r a ds i(- U -(-a+l)i)logP(T i _i,T i ) 



a 2 — a — v 2 + i (—2a + 1) v 



(52) 



(53) 



In the denominator of the expression above, we recognize the time t conditional discounted 
characteristic function of the zero coupon bond at time T,_x having maturity Tj evaluated at the 
point v — (—a + l)i. While this object may seem difficult to compute, we remarkt that under our 
assumption on the shape of the short interest rate 

r t = h d +Tr[H d j:(t)] , (54) 

we can perform the computation, by relying on the techniques we have been already employing. 

Proposition 8. The time t conditional discounted characteristic function of the logarithmic zero- 
coupon bond price with maturity Ti, at time Ti—i, evaluated at the point U3 is given by 



e - Sl^ 1 {h d +Tr[H d Y,{S)])ds+u log P(T^ 1 ,T i ) 

= exp {u>A zc {Ti - Ti_i) + A LZC (Ti_i - t) 
+Tr [A LZC (T i .. l -t)Z{tj\} 

where the deterministic functions A LZC ,B LZC satisfy the system of matrix ODE 

Q A LZC 



dr 
8B LZC 
dr 



Tr [/3Q t QB lzc (t)] - h d , A LZC (0) = 0, 



= B LZC (r)M^ + 1 B^(t) 



T nLZC i 



+ 2B lzc (t)Q t QB lzc (t) - H d , B LZC (0) = U)B zc (Ti - Ti-i), 
whose solution is give by 



A lzu (t) 



-Tr [log B^ zc (r) + rM^]-h d r, 



(55) 



(56) 



(57) 



(58) 



B lzc (t) = {uB zc {T % T^B^i^+B^ir))- 1 



X [L}> 



B zc {T t - Ti-{)B^ c {j) + B^ zc (t)) 



and B% z (t), B% z {t) are submatrices in: 



BtfC{r) B^ zc {r) J- 6XP 



-2Q T Q 



-H d -M Qd ' T 

With the previous result, caps and floors can be easily priced within our framework. 



(59) 



(60) 



6 The multi-Heston model of De Col et al. (2011) as a 
nested approach 



In this section we follow [Benabid et al. (20101 in order to show that when all the matrices 
M, Q, R, S(0), Ai, H l are diagonal and for a particular specification of the noise Z driving the 
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exchange rates, the Wishart model becomes equivalent to the multi-Heston approach introduced 
in De Col et al. (2011 ). In this particular case the dynamics of the elements of the Wishart process 
take the simpler form 



and 



cE pp (t) = (pQ 2 pp + 2M ii Z pp (t)) dt + 2Q PP V \/W) pk dWg{t) 



fe=i 



We notice that 



k=l 



dZ pm (t) =S pro (t) (M pp + M mm ) dt + Q pp J2 VW) km dW kp (*) 

d 



fe=i 



so the d-dimensional process W = ( W(t)} , with 

v / t>o 



\t) = ^ p - p J2VW) pk dwg 



kp(t) 



k=l 



is a vector of independent Brownian motions. As a consequence, we may conveniently express the 
dynamics of a generic diagonal element as 

dH pp (t) = {(3Q 2 pp + 2MuSp P (t)) dt + 2Q pp ^ pp (t)dWf (t). 

Assume also that the matrix Z^ — (^Z® (t)j is diagonaj^j A quick inspection at the 
diffusion part of the dynamics of the exchange rate S l,J reveals that we are considering a setting 



that is equivalent to the multi-Heston model of De Col et al. (2011 ), namely, we may conveniently 
rewrite the dynamics 



as 



dS^(t) 
6W(t) 



dS^jt) 
S^(t) 



= (r l (t) - r j (t))dt + Tr \{A, - Aj)y/E(fjdZ Q ' \t) 



(r l (t) - r ] {t))dt + (a, - aj -) T yj Diag (Y,(t))dZ q ' (t), 



where a^a^ are column vectors, Z^ (t) is a d-dimensional vector Brownian motion with Z 1 ^ (t) = 

Z®p(t) and \jDiag (E(i)) denotes a diagonal matrix, featuring the elements y/T, pp (t), p = 1, d 
along the main diagonal. The short rate processes may be equivalcntly expressed by means of a 
scalar product between a vector featuring the elements of the main diagonal of H l and a second 
vector featuring the elements of the main diagonal of £. This allows us to connect the finite 
dimensional distributions generated by the hybrid Wishart process with those generated by a 
multi-Heston model. Being the finite dimensional distributions equal, also the prices of European 
options will be the same under the two specifications. From this informal discussion, we also 



deduce that the Wishart hybrid model nests an extension of the approach of De Col et al. (2011) 
featuring stochastic interest rates. 



Note that Z is no more a matrix Brownian motion, which by definition consists in d X d independent Brownian 
motions. 
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7 Numerical illustration 



7.1 Joint calibration of the FX-IR hybrid model 

In this subsection we report the results of a simultaneous calibration of a foreign exchange options 
volatility surface and of the two yield curves of the economies linked by the foreign exchange rate. 
To be more specific, we consider an implied volatility surface of EUR/USD and the yield curves 
of EUR and USD economies. 



7.1.1 Description of the data 

We consider market data on January 15 th 2013. The data set features an implied volatility surface 
for options written on EUR/USD. Moneyness ranges from 5A put up to 5A call. Maturities 
range from 1 day till 15 years. The conversion between deltas and strikes can be easily performed 
along the lines of Clark (2011 ) or Wystup and Reiswich (20101. As far as the interest rate market 
is concerned, we obtained the yield curves for both EUR and USD. In this case we considered 
maturities up to 20 years for both curves. 



7.1.2 Calibration results 

The idea of the present calibration is to try to fit simultaneously the following data: the market 
FX implied volatility surface of EUR/USD and the EUR and USD yield curves. Our calibra- 
tion is thought of as instrumental for the evaluation of long-dated FX products like power reverse 
dual currency notes (PRDC, see the introduction). For this reason, the range of expiries that we 
consider is much larger. As far as the yield curves are concerned, we fit market data up to 20 
years. Concerning the implied volatility surface, we consider maturities ranging from 1 month up 
to 15 years (the longest maturity at which traders are allowed to trade). The penalty function is 
constructed first by considering the distance between market and model implied volatilities for FX 
options and then by also looking at the distance between market yields and model yields computed 
according to formula (|49|). 



Figures [T] and [2] report the results of our calibration. The fit is very satisfactory across differ- 
ent strikes and maturities for FX options. The upper part or Figure [T] shows the actual distance 
between the two surfaces, that are almost overlapping. For the sake of readability in the bottom 
left part of Figure [T] we multiply by a factor > 1 all model implied volatilities so as to ease the 
comparison. The bottom right part of Figure [l] provides a more precise view on the quality of the 
calibration of the FX surface, by plotting for each point in the maturity/delta space the squared 
difference between model and market implied volatilities. Figure [2] reports a comparison between 
market and model implied yield curves. Red stars and blue circles denote model and market im- 
plied yields respectively. Even in this market we are able to obtain a satisfactory fit. Recall that 
we perform a joint calibration procedure, meaning that the same set of model parameters allows 
to obtain the results of Figures [l] and [2] 

The values of the parameters arising from the calibration are reported in Table [T] By observing 
the values of the parameters we obtained, we notice that for the parameter j3 we have /? > d + 1, 
hence for the fitted model we conclude that the dollar and the euro risk-neutral measures are 
well posed martingale measures. As far as the values of Heur and hjjSD are concerned, they 
are negative, which is not in line with our starting assumption. In a first calibration experiment, 
we imposed the constraint hi > 0, i = USD, EUR but observed that the model was not able 
to replicate the observed yield curve shape. Recall that in the present framework we are using 
(projections of) the Wishart process for the fit on both the FX implied volatility surface and the 
yield curves. This fact results in a trade-off between the parameters of the model. In other words, 
a small value of the initial state variables E conflicts with the ability of the model to fit the short 



term smile. This is in line with the findings of Chiarella et al. (2010) who also found negative 



values for h, thus implying a distribution of the short rate that can become negative. Anyway, 
upon relaxation of the positivity constraint we obtained the very satisfactory fit that we reported 
in the present section. 
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8 Conclusions 



In this paper we introduced a novel hybrid model that allows a joint evaluation of interest rate 
products and FX derivatives. The model is based on the Wishart process and so it retains the same 
level of analytical tractability typical of the affine class, like in the multi-Heston version introduced 
in De Col et al. (2011 ). In analogy with De Col et al. (2011 1, we have a model that is theoretically 



consistent with the triangular relationship among FX rates and other stylized symmetries that 
are commonly observed in the FX market. Moreover, the model allows for interest rate risk since 
there is the presence of stochastic interest rates that are non trivially correlated with the volatility 
structure of the exchange rates. A successful calibration of the EUR/USD implied volatility sur- 
face and of the EU R and USD yield-curves suggests that the framework we propose is a suitable 
instrument for the evaluation of long-dated FX products, where a joint description of interest rates 
and FX markets is required. 

There is ample room for future research. For example, the introduction of Sf valued jump 



processes, as in e.g. Muhle-Karbe et al. (2012), would help in capturing the highly skewed implied 
volatilities that are usually observed for short maturities. Another interesting direction is the 
generalization to the multi-curve setting that is emerging in interest rate modelling as a response 
to the recent financial crisis. Several authors have attempted a solution in the FX context in order 
to model different levels of risks between e.g. EONIA and EURIBOR curves (see e.g 



( |2010| , |Fries| (|2010j), |Kenyon| ( |2010[ ), |Kijima et al.| (|2009j) and references therein). We believe that 
a Wishart based approach would give very interesting insights in this multivariate puzzle. 



Bianchctti 



9 Proofs 

9.1 Proof of Proposition [T] 

Application of the Ito formula to the product S l ' L (t)S L > : > (i) (using the property Tr[dW^4]Tr[dH^B] 



dS^ l (t) 



Tr[AB]) leads to Q. In formulas: 

= ( r < _ r l )dt + Tr[(A z - Ai)V{t)Ai]dt + Tr[(A t - A{)^Z(f)dZ{t% 



S*.'(t) 
d 9' J (t ) 

-^TJJjy = {r l -r^dt + TrliM-A^Mdt + TrUt-A^^n^dZit)], 

dS^{t) =dS l > l {t) t S\> 3 + S l \t) t dS\> 3 +d(S i ' l ,S l ' j ) t 

=Sl> j S i > l (t)((r i - r l )dt + Tr[(Ai - 4)E(f) A^dt + Tr[(Ai - A^^tjdZit)} 



+ S l t 3 S l ' l (t) (r l - r 3 )dt + Tr[{Ai - Aj)H{t)Ai\dt + Tr[{Ai - Arfy/ftydZit)] 



d( I Tr[(A l -A l )^/z(u)dZ(u)}, Tr[{Ai - A 3 )^{u)dZ{u)} 
lo Jo It 
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We concentrate on the covariation term: 



d( / Tr[{Ai - A l ) y / JMvjdZ(v)}, / Tr[(A, - Aj)y/nv)dZ(v)\ 
wo JO / 

/ p. d p. d 

=d( Y (Ai-A^^yj^dZiv)^, / Y (Ai-A/) st ^) tu dZ(v) us 



o,q,r— 1 



p,q,r,s,t,u=l 

x (Ai - A 3 ) st y / T,(t) tu dZ(t) us 5 r=u S p=s 

d 

- Y ( A i- A iUV^) q ui A i- A j)stVW) tu dt 

s,q,u,t—l 
d 

= y ( A i- A iUVm q uVm u t( A i - A J) ts dt. 

s,g,u,i— 1 

By assuming that the matrices Ai,Ai,Aj are symmetric we get 



Tr[{A t - AfijnyjdZiv)}, / Tr[(At - A j )y/Z(y)dZ(v)] 

Jo 

= Tr[(A i -A l )E(t)(A,-A J )]. 
Finally, using the fact that terms in the trace commute we obtain: 

= {r i _ r j )dt + Tr[{A t - A^Mdt 

+ Tr[{A t - Aj)y/nftdZ{t)} 



as desired. 



9.2 Proof of Proposition [4] 

Recall that <\> % >■? (w, t, r, x, E) = G IJ («J, i, r, x, E) where G was previously defined as G 1 '^ (ui, t, t, x, E) 
gQ jg^^Tj^ x l 'j{t) := log S t 'i(t). These functions represent the conditional discounted character- 
istic and moment generating functions of the log-exchange rate respectively. We follow closely 



Da Fonseca et al. ( 2008 ) in order to determine these quantities, so we first write the PDE satisfied 



by G, which requires the dynamics of x(t) = x 1,3 (t) under the measure Q l : 

d\ogS^(t) = (V - r j ) - l -Tr [(A t - Aj) E(i) (A t - Aj)]j dt 



+ Tr 



(Ai - Aj) y/E(t)dZ Qi (t) 



where the short rates are driven by the Wishart process in line with 
solves the following PDE in terms of r = T — t: 

d 



G IJ = Ax S G !J - r l G l 



G 4J (o;,r,0,a;,E) = e" 



(61) 

The Laplace transform 

(62) 
(63) 



To solve this PDE we first determine the infinitesimal generator A x ,s ■ This will feature the contri- 
bution of three terms: the process x, the Wishart process E and the mixed term that corresponds 

pi 2 

to the coefficient of the term a 1^ — and arises from the correlation structure. The first is trivial, 
the second is known thanks to Bru ( 1991 ) and is of the form: 

Tr 



(fin T + ME + EM T j D + 2Y,DQ T QD 
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where D is the differential operator: 



D 







<9£ 



pi 



In order to compute the mixed term, we notice that under the measure 
d(logS i > j ,X pt ) u 



M \k Tr 



q,r,s—l 



2d( [' Yl ( A i- A j) ab VW) bc dZZ(v), [ J2 VE pq (v)dZ$(v)Rj s Q st ) S c=q S a= 

\J0 Jo / 



q y r,s=l 



= 2 ]T (A i -A j ) rb ^/m bq Vm pq RlQstdt 

b,g,r,s— 1 
d 

= 2 J2 Z(t) P b{A~A 3 ) br Rj s Q st dt, 

b,q,r,s—l 

where we have used the fact that 

2Tr [E (Ai - Aj) R T QD] — = 2 £ A P S p6 - A,-) fer i?^ s Q st 



<9x 

p,b 7 r,s— 1 

and that D is symmetric. Now we can state the PDE satisfied by G: 
dG 



dx 



j ) - ^ Tr - Aj) £ (Ai - Aj)] \ ^ 



l T r[(A i -A j )E(A i -A j )}^ 



Tr 



nn T + ME + EA/ T ) DG + 2 (ZDQ T QD) G 



dG 



+ 2Tr [S - Aj) i? T QL>] h l - Tr [iPE] . 

Since the Wishart process £ is affine we guess the following 

G 1 - 3 (w, t, t, x, S) = exp [C(r)x + _4(r) + Tr [B(r)S]] , 



(64) 



(65) 



with C,A EM. and B € Sd s.t. the transform is well defined, moreover these functions satisfy the 
following terminal conditions: 

A(0) = 0el. 
C(0) = w€M, 
B(0) = 0eS d . 



We substitute the candidate ( 65 ) into ( 64 ) and obtain 

d 



Or 



C(t) = 0, 



hence: C(t) = cj Vr. We also have the following (matrix) Riccati ODE: 
d 



8t 



B = B(t) ( M + loQ R (Ai — Aj) ) + (M +u (A, -A j )R , Q) B(t 



2B(t)Q t QB(t) 



(Ai - Aj) 2 + (w - l)H l ~ uH j 



(66) 
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and the final ODE which may then be solved upon direct integration 

a 



A = oj (h l - h j ) - h l + Tr [nn T B(r)] . 



Following Grasselli and Tebaldi ( 2008| ) , it is possible to linearize ( 66 ) by writing 

for F(t) € GL(d) and Q(t) £ Md, then we have 
d_ 

d_ 



F=-F{t)[AV +Lj(A l ~A j )R , Q) -2G(t)Q'Q 



Q = G(t) [M + ujQ R (Ai - Aj 



^—^ {Ai - Aj) 2 + (« - l)H l - wH* 



with J-(0) — Id and G(0) = 6(0). The solution of the above system is 



x expr 



(B(0),I d ) 



Bii(r) B 12 (t 
B 21 (t) B 22 (r 

M + ioQ T R(Ai - Aj) 



-2Q T Q 



(A, - Aj) + (u- i)H l - loW - { M 1 + lj (Ai - Aj) R 1 Q 



so that the solution for B(t) is 

B(t) = (B(0)B 12 (t) + B 22 {t))- 1 (B(O)Bu(r) + B 21 (r)) . 
Since B(0) =0we finally get B(t) = B 22 (t)~ 1 B 2 \(t). Now starting from ( [69] ) we write 
1 / d 



2\dr 



T + T{t) (m t +u(Ai-Aj)R T QU (Q T Q) 1 =Q(t) 



We plug this into (68), then we insert the resulting formula for B(t) into (67 1 and obtain 



a 

dr 



A = lu (W - h j ) - h* 



Tr 



^^(t)- T + (M T + u (Ai - Aj) R T Q 



whose solution is 



A=(u (h* - h j ) - fj) t - ^Tr 



log J-(r) + (M T +cj(Ai — Aj) R T Q 



9.3 Proof of Proposition [5] 



(67) 



(68) 



(69) 



(70) 



The matrix Riccati ODE ( 66 ) may be rewritten as follows after replacing Q by a small perturbation 
aQ, a £ K: 



d_ 



B =B(t) (m + aujQ T R (A t - Aj)) + (M T + uoa [A, - Aj) R T Q S j B{t) 



2a 2 B(r)Q T QB(t) + (A t - Aj) 2 



B(0) =0. 



(71) 
(72) 
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We consider now an expansion in terms of a of the form B = B° + aB 1 + a 2 B 2 . We substitute this 
expansion and identify terms by powers of a. We obtain the following ODE's. 



(73) 



dr 


=B°(r)M4 


M^(r) + W2 - W (^ Aj) 2 


—B 1 
dr 


=6 1 (t)M4 


M T B 1 (r) + B°{t)Q t R (Ai -Aj)u 


U 2 
dr 


+ uj(A 2 - 


Aj) R t QB°(t) 


=S 2 (t)M4 


M t B 2 (t)+B 1 (t)Q j R(Ai -Aj)u 




+ oj(A 1 - 


A,) R t QB 1 (t) + 2B a (r)Q T QB a (r) 



Let 7 



then these equations admit the following solutions 



B°( r ) / e (T - u)M (Ai - A,) 2 e^ M du 

2 Jo 

:=7B°(r), 

B x (t) =tw jT r e( T -^* T (tf»Q T i?(A 4 - Aj) 

+ (Ai - Aj) i? T Q6°(u)) e^ T ^ £l du 
:=7wB 1 (r), 

6 2 (r) = 7 ^ 2 jT e (-«)^ T (b 1 {u)Q t R (Ai - Aj) 
+ (A, - Aj) i? T Q6 1 (u)) e^- u ^ [ du 



+ 7 2 / e (-«) M ' 2S°( u )g T Q6°( W ) e ( r - u)M d W 
Jo 

= 7 2 S 20 (t)+ 7W 2 6 21 (t). 



(74) 
(75) 



(76) 



(77) 



(78) 



whereby we implicitly defined the matrices B° (t) , B 1 (t) , B 20 (t) , B 21 (t) . We can now write the 
function B(t) as follows 



B(t) = 7#°(t) + a 7 uB 1 (r) + oVB 20 (t) + a 2 7 w 2 £? 21 (r). 



(79) 



A direct substitution of ( 79 1 into <j67J) allows us to express the function A(t) as 

nn T [ B°(u)dv 
Jo 



A(t) =lu (r, — rj) t + ^Tr 



ajujTr 



nn T [ B x {u)dv 
Jo 



a 2 7 2 Tr 


nn T [ B 20 (u)du 


+ a 2 joj 2 Tr 


nn T f B 21 (u)du 




Jo 




Jo 



=u (n -rj)r + 7A°(r) + orywi 1 ^) 
+ a 2 7 2 A 20 (r) + aV 2 A 21 (r), 



(80) 



having again implicitly defined the functions A°(t), A 1 (t), A 20 (t), A 21 (t). We consider now the 
pricing in terms of the Fourier transform, i.e. uj = iX, as in (24). Let Z denote the strip of 
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regularity of the payoff. A Taylor-McLaurin expansion w.r.t. a gives the following: 



C(S(t),K,r) 



2tt 
fa 



i\(r t -r 3 )T+i\x+ 1 (A a {T)+Tr[B°(r)^})^^ dx 



2tt 



(A 1 (r) + Tr [B 1 ^ 

7 iAe iA(r '^ )r+iA2;+7 (^ 0(T)+Tr ^ 0(T)E ])$(A)dA 



(A 2Q {T) + Tr [s 20 (r)S 

— / 7 2 e iX{r *- r i )T+iXx+ ^( A °( T)+Tr i® 0{T) ^)<I>(\)d\ 



a 



2tt 

2 



(a 2 \t) 



Tr 



7iA 2 e iA(n-r 3 )r+iA, ;+7 (j; (r)+Tr[e (r)S]) $( - A ^ A 



(i 1 (r)+Tr [e 1 (r)S 



2tt 



y 2 i 2 A 2 e iA(r I -r 3 )r+iAx+ 7 (^ (r) + Tr[g°(r)S]) $ ^ A ^ A ^ 



Recall now from (40) the definition of the integrated Black-Scholes variance. In the previous 
formula in the first term we recognise the Black Scholes price in terms of the characteristic function 
when the integrated variance is v = <j 2 t: 



C B &s (S(t),K,a,r) 



2tt 



iA(r^ — Tj )T J ri\x-\- — 



(81) 



so that the price expansion is of the form: 

C(S(t),K,r) ^C BSz s(S(t),K,a,T) 

+ a (A 1 (r) + Tr W (r)sl ) d 2 xv C Bhs (S(t), K, a, r) 



+ a 2 (A 20 (T) + Tr £ 20 (r)£ ) d 2 2 C BSiS (S(t), K, a, r) 
-Tr 

Tr 



a 2 (A 21 (t) 



B\t)t\) 2 d^Cs^s (S{t),K, a, t) . 



(82) 



which completes the proof. 



9.4 Proof of Proposition [6] 

We follow the procedure in |Da Fonseca and Grasselli (20111. We suppose an expansion for the 
integrated implied variance of the form v = <J 2 mp T = Co + o-Ci + and we consider the Black 
Scholes formula as a function of the integrated implied variance and the log exchange rate x = 
log 5: CB&s(S(t),K,&, T ) = CB^s(x{t),K,a 2 mp T,T). A Taylor-McLaurin expansion gives us the 
following: 

CB&s(x(t),K,crj t,t) = C B & s (x(t),K, Co, t) + a(id v C B &s(x(t),K, ( ,t) 



a 



+ — (2{ 2 d v C Bg6S (x(t),K,Co,T) 
+(ld 2 v2 C B s.s(x(t),K 1 (o,T)). 



(83) 
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By comparing this with the price expansion ( 82 ) we deduce that the coefficients must be of the 
form: 



Co = v Q 

Ci 

C2 





B x (r)S 


j dl v C B hs 







(84) 
(85) 



-Ci 2 df 2 C B&s + 2 (A 20 (t) + Tr [s 20 (t)e]) d 2 v2 C B&cS 



(A 21 (t) 



2d v C B &s 
B 21 {T)Y\)dl 2v C Bhs 



Tr 
2d v C B iks 



A 1 (r)+Tr S^t)^ 



d\i,.iC B us 



2d v C B us 



(86) 



where the Black Scholes formula C B &s(x(t), K, of mp T, r) is evaluated at the point {x, K, vq,t). In 



order to find the values of Ci>C2j we differentiate ( 32 )-( 35 ) thus obtaining the following ODE's: 
d 



B v = S u (r)M + M T B°{r) + (A, - Ajf , 



ch 



■^B 1 =B\t)M + M t B\t) 



B°(t)Q 1 R (Ai - A,) + {A, - Aj) R 1 QB v (t) 



d ~ 



B M = B m {t)M + M t B 20 (t) + 2B v {t)Q 1 QB u {r), 



Ot 



^B 21 =B 21 (t)M + M t B 21 (t) 

+ B\t)Q t R (Ai — Aj) + {At - Aj) R t QB\t). 
We consider a Taylor-McLaurin expansion in terms of r 

B°(r) « (A, - Aj) 2 r + ^ ((A - A 3 ) 2 M + M T (A, Aj)' 

2 

&(.t) « ^- [(A* - Aj) 2 g T i? (A -Aj) + (A - A,-) i? T Q (A - A,- 



2 

l! 
6 

M 1 

T l 
6 



(A - Aj) 2 Q T R (A - Aj) + (A - A,) R 1 Q (A - A,)' 
(A - A,) 2 Q T i? (A - A,) + (A - Aj) R t Q (A - A,) 2 ] ) 



M 



(A - AjY M + M 1 (A - A,-) 2 Q T i? (A - A 3 ) 



+ (A - Aj) R T Q [(Aj - A,) 2 M + M T (A, - A,-)' 

6 20 (r) « ^4 (A ~ A,) 2 Q T Q (A - A 2 f 

B 21 (r) » ~ ( [(A ~ A,) 2 Q T i? (A - A,) + (Ai - A,) R T Q (A, - A,-) 5 
x Q T R (A, - A 3 ) + (A, - Aj) i? T Qx 
(A - A,) 2 Q T i? (A - A 3 ) + (A, - Aj) i? T Q (A - Aj) 2 



(87) 



(89) 



(90) 



(91) 



(92) 
(93) 



(94) 



Noticing from (36)- (39) that A 1 (t) are one order in t higher than the corresponding B 1 (t), the 
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following approximations hold: 

A (r)+Tr [b°(t)S 
A\r)+Tr Ib 1 ^ 
A 20 (r)+Tr fe 20 (r)S 
A 21 (r)+Tr fe 21 (r)E 



= Tr 



(Ai - Aj) E 



o(r) 



Tr 

2 
3 
Tr 



(A i -A J ) 2 Q T i?(A l - Aj-)S t 2 +o(t 2 ) 



= ^Tr 



(A-^) 2 Q T i?(A-A,) 



+ o(r 3 ) 



(95) 
(96) 
(97) 



(Ai- Aj)R T Q(Ai- Aj 



(98) 



and the variance V 



x Q T R(A t - Aj-) eJ — + o(t 3 ). 

We introduce two variables: the log-moneyness m/ = log ^- — — 

Tr [(Ai — Aj) E(t) (A, — Aj)] r. Then, from Lewis (20001, we consider the following ratios among 
the derivatives of the Black-Scholes formula: 

9l v C B SzS (x,K,V,t) 



d v C B &,s (x, K, V, t) 
d 2 v2 C BglS (x,K,V,t) 
d v C B k.s (x, K, V, t) 
d 3 x 2 iV C BSlS (x,K,V,t) 

d v C B gzS (x,K,V,t) 
dj^CmzS (x,K,V,t) 
d v C B &s (x,K,V,t) 



1 



1 m/ 

2 V 

m ) 

2V 2 2V 

1 rrif — 1 
4 + V 
rrif m 2 (m,f — 1) 
W 1 + 2l^ ' 



V 2 ' 



(99) 
(100) 
(101) 
(102) 



Upon substitution of (|95|)-(|102|) into (|85|)-(|86|) we obtain the values for i = 0, 1, 2 allowing us to 

(103) 



express the expansion of the implied volatility. 

Co = Tr [(Ai - A 3 ) E(i) (Ai-Aj)]r, 

(A l -A J ) 2 Q T R(A l -A J )^(t) 



Tr 



Ci 



Tr[(A-^)E(i) (Ai - Aj)} 

,2 

rT 



(104) 



7 



Tr [(A - Aj) E(t) (A - A^ 



1. 



Tr 



(A l -A,) 2 Q T Q(A-A J ) 2 E(t) 



+ -Tr 



(A - Aj) 2 Q t R (A - Aj) + (A l - Aj)R> Q (A, - Aj) 



x Q T R(Ai - Aj)T,(t) 



Tr 



(A- A j ) 2 Q T R{A i - Aj)Z(t) 



4 Tr [(Ai - Aj) S(i) (A; - Aj)} 
By plugging these expressions we obtain the result. 



(105) 
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10 Images and Tables 



Parameter 


Value 


Parameter 


Value 


£(t)(l,l) 


0.1688 


A us (l,l) 


0.7764 


£(t)(l,2) 


0.1708 


A us (l,2) 


0.4837 


£(t)(2,2) 


0.3169 


^4ms(2, 2) 


0.9639 


M(l,l) 


-0.5213 


^4e«r(lj 1) 


0.6679 


M(l,2) 


-0.3382 


A eur (\, 2) 


0.6277 


M (2, 1) 


-0.4940 


A eU r (2, 2) 


0.8520 


M(2,2) 


-0.4389 


hus 


-0.2218 


Q(l,l) 


0.2184 


heur 


-0.1862 


Q(l,2) 


0.0957 


H us (l,l) 


0.2725 


Q(2,l) 


0.2483 


H us (l,2) 


0.0804 


Q(2,2) 


0.3681 


H us (2,2) 


0.4726 


J2(l,l) 


-0.5417 


Heur ( 1 ) 1 ) 


0.1841 


i2(l,2) 


0.1899 


H eur {l,2) 


0.0155 


i?(2,l) 


-0.1170 




0.4761 


R(2,2) 


-0.4834 


P 


3.1442 



Table 1: This table reports the results of the calibration of the FX-IR Wishart hybrid model 



27 



USDEUR - Calibration result 




Figure 1: We consider a joint calibration of the implied volatility surface of USDEUR and of the 
two yield curves of USD and EUR simultaneously. The figure on the top compares the implied 
volatility surface to observed market data. As the fit is very good and the two surfaces are 
almost indistinguishable, on the bottom left figure we shift upwards the model implied volatility 
by multiplying each point by 1.1 so as to encrease readability. On the bottom right we report 
squared errors in implied volatility for each point of the surface. 
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